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The Rayleigh diffraction bound sets the minimum separation for two point objects to be distin- 
guishable in a conventional imaging system. We demonstrate resolution enhancement beyond the 
Rayleigh bound using random scanning of a highly-focused beam and A^-photon photodetection 
implemented with a single-photon avalanche detector array. Experimental results show resolution 
improvement by a factor ~\/N beyond the Rayleigh bound, in good agreement with theory. 
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The response of a diffraction-limited imaging system 
to a point-like source — its point spread function (PSF) — 
has extent inversely proportional to the entrance pupil's 
area. The image is obtained by convolving the system 
PSF with the light distribution of the object. There- 
fore, details finer than the PSF's extent are lost under 
conventional (entire-object) illumination. The Rayleigh 
diffraction bound sets the minimum separation for two 
object points to be distinguishable in the image. 

Two classes of quantum strategies have been suggested 
to circumvent this bound. The first relies on techniques 
from quantum metrology pj] , in which the image informa- 
tion is encoded into suitably tailored non-classical light 
beams [^tJI- For example, quantum lithography [5| ex- 
ploits the effective de Broglie wave-function of TV pho- 
tons in a delicately crafted state to obtain an increase in 
resolution proportional to N (see [6| for recent develop- 
ments). These methods are typically extremely sensitive 
to photon loss and noise, because they rely on delicate 
quantum effects (such as squeezing and entanglement). 
Thus they are best suited to short-distance applications, 
such as microscopy, where losses can be controlled, as op- 
posed to standoff sensing, such as laser radar operation 
over km or longer path lengths, for which substantial 
diffraction and atmospheric losses will be present. 

The second class of quantum strategies for beating 
the Rayleigh diffraction bound exploit postselection [7J 
to extract the high-resolution image associated with a 
non-classical component from classical-state light con- 
taining information about the object to be imaged ISl- 
17| . Because postselection involves discarding part of 
the measurement data, these procedures inherently suf- 
fer detection inefficiency that increases the time required 
to acquire an image. However, their spatial resolution 
can nonetheless exceed the Rayleigh diffraction bound. 
Furthermore, because they employ classical-state (laser) 
light, these techniques degrade gracefully with increasing 
loss and noise, making them suitable for standoff sensing. 
In this Letter we report the first experimental demonstra- 
tion of one such technique, viz., that of Giovannetti et al. 
[17j , in which the object is illuminated by a focused light 
source and scanned. The scanning pattern is irrelevant 



so long as the area of interest is covered, i.e., a random 
scan pattern suffices. The image is formed using only 
those pixels that count exactly N photons within a mea- 
surement time T. The expected resolution inrprovement 
is \/iV over standard entire-object illumination, until the 
limit set by the focused-beam illumination. We begin 
with the theory for this technique. 

Theory: — We are interested in an active imager, such 
as a laser radar, comprised of a transmitter and a receiver 
in which we control the object illumination and form the 
image with the receiver. For such systems, the spatial 
resolution is a function of two antenna patterns, viz., the 
transmitter's illumination pattern on the object and the 
receiver antenna pattern, set by the diffraction limit of its 
optics, projected onto the object. When floodlight illumi- 
nation is employed, so that the entire object is bathed in 
light, the resolution limit is set by the receiver's Rayleigh 
diffraction bound. If the transmitter and receiver are 
co-located, they can share the same optics so that their 
antenna patterns have identical Rayleigh bounds whose 
product gives the overall resolution behavior. Alterna- 
tively, if a small transmitter is located much closer to the 
object than is the receiver, it is possible to project very 
small spots onto the object to be imaged. If this is done 
in a precision scan, so that the receiver knows exactly 
where the transmitter is pointing at any instant in time, 
a simple energy measurement at the receiver will realize 
resolution limited by the transmitter's antenna pattern, 
regardless of the receiver's own Rayleigh bound. How- 
ever, creating that precision scan, and relaying the scan 
positions to the receiver, could easily be a major chal- 
lenge, especially if the transmitter is mounted on a small 
unmanned air vehicle. The postselection technique of 
Giovannetti et al. circumvents that problem by allowing 
the scan pattern to be arbitrary, even random, so long as 
the object's region of interest is covered. 

Suppose that a focused transmitter emits a +z- 
going, quasimonochromatic, paraxial, linearly-polarized 
laser pulse with scalar complex envelope ET{p,t;6) — 
.yiNT/nD^ g(^)g-«fe|p|V2LT+»fce-p fo^ \p\ < Dt/2, where 
p = {x, y) is the transverse coordinate vector, 6 = 
{(^x,Oy) is the transmitter aim angle, and k is the wave 



number at the center wavelength A. We will normal- 
ize Et to have units ^photons/m^s, and take the pulse 
shape s{t) to satisfy Jdt |s(i)P = 1, so that Nt is the av- 
erage transmitted photon number |18j . This pulse tran- 
silluminates an object Lt-'ch away from the transmitter 
|19{ . The light that passes through the object is then 
collected by a diffraction-limited circular lens of diame- 
ter Da located L/j-m in front of the object. The focal 
length of this lens is such that it casts an image of the ob- 
ject at a distance Lj > Lr beyond the lens. Neglecting 
the propagation delay and correcting for image inversion, 
the photon-flux density in the image plane is then 

|i?iM(PiM,i;e)|'-iVT|s(i)|'x 
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where 0{p) is the object's field-transmission function, 
m = Lj/Ln is the image magnification 20], and we see 
the transmitter and receiver's circular-pupil antenna pat- 
terns. 

To exhibit the sub-Rayleigh resolution capability of the 
scheme from [17| we shall assume that: (1) the transmit- 
ter's antenna pattern fully resolves all significant features 
in 0{p) [2l|; (2) the image-plane photon counting array 
has pixels of area Ap sufficiently small that diffraction, 
rather than pixel size, limits image resolution; and (3) 
the photon counting array outputs pixel counts taken 
over the full T-s-long extent of s{t). For a given illu- 
mination angle 6, the pixel counts are then statistically 
independent, Poisson random variables with mean value, 
for the pixel centered at pj^, given by 

Ji{TrDB.\9LT - Pim/^I/^Lr] 
nDR\9LT - PIM/"^|/2Aii^ 

where rj is the detector's quantum efficiency. If photon 
counts are collected from the pixel at pi^ while 6 is ran- 
domly scanned over the object region, the unconditional 
probability of getting N counts from that pixel is 
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where p{6) is the scan pattern's probability density func- 
tion. Postselecting those pixels for which N counts have 
been registered, we get an image Zjv(pji^) ex -P/v(Pim)- 
For N > maxNg{pij^), the Poisson distribution is mono- 
tonically decreasing with increasing N, whence 
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where we have suppressed multiplicative constants and 
ignored the exponential term as it is independent of N . 
Here we see that the postselected image contains jOp^ 
convolved with the iVth power of the receiver's Airy disk 
pattern, i.e., a point-spread function that is ^yN nar- 
rower than the Airy disk itself when compared on a main- 
lobe area basis |17|. Note that violating the preceding 
monotonicity condition can lead to the "donut-hole" ef- 
fect exhibited in the experiments described below. 

Experiment: — We demonstrate the concept of sub- 
Rayleigh imaging with the setup shown in Fig.[lja). The 
object to be imaged in transmission was part of a U.S. 
Air Force (USAF) resolution target consisting of alter- 
nate opaque and clear stripes of width 125 /im (4 line 
pairs/mm), as indicated by the arrow in Fig.[2ja). A 532- 
nm laser was mounted on an XY translation stage that 
provided scan coverage over the entire object with a 20- 
/im-radius focused spot. We imaged the object through a 
/ = 25-cm diffraction-limited lens set in a 2-mm-diameter 
aperture. The optics provided 5.3 x image magnification, 
yielding a 660-/xm-wide strip at the image plane. Under 
conventional (entire-object) illumination, shown in the 
setup of Fig. [Ijb) , the Rayleigh diffraction bound for the 
imaging system at the image plane was 1.86 mm, which 
is ~2.8x larger than the stripe width. Figure[2][b) shows 
the conventional-illumination image that was obtained 
using standard photodetection (with all events counted) : 
the stripes are unresolved, as expected. Note that with 
full illumination of the object, we were not able to go 
beyond the Rayleigh bound even with the A^-photon de- 
tection scheme, indicating that focused illumination is a 
necessary requirement. 
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FIG. 1: (color online). Setup schematics for (a) sub-Rayleigh 
imaging with focused illumination, and (b) conventional co- 
herent imaging with full illumination. 

The detector was a compact 32 x 32 Si single-photon 
avalanche diode (SPAD) array fabricated with a com- 
plementary metal oxide semiconductor (CMOS) process 
[22| . The Si SPAD was a p-n junction reverse-biased 



above its breakdown voltage and operated in the Geiger 
mode with a detection efficiency of ~30% at 532 nm. 
Each pixel of the CMOS SPAD array consisted of one 
SPAD with its front-end active quenching and resetting 
electronics and a digital counting circuitry for in-pixel 
pre-processing. The pixel pitch for the array was 100 /im, 
and the SPAD had a fill factor of only 3.1% at each pixel 
due to the presence of on-chip electronics. Owing to the 
large separation between SPADs we did not observe any 
cross talk. The average pixel dead time (including all the 
electronic circuitry) was 300 ns. Each SPAD delivered 
a digital output pulse for every single-photon detection 
event with no readout noise. The in-pixel counting cir- 
cuitry would compute the number of single-photon events 
within its user-selectable integration time of 1 /is or more 
and store the tally in an in-pixel memory cell. By mea- 
suring incident photons over a long integration time, N- 
photon sensitivity in the time domain can be achieved at 
the single-pixel level. The array readout was performed 
through an 8-bit data bus without interrupting the next 
1024-pixel frame of photon-counting integration, and the 
maximum frame rate was 10^/s. A typical integration 
time of each frame was tens of fis. 
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FIG. 2: (color online), (a) Object indicated by arrow: tar- 
get # 1, group 2 of a USAF resolution target, composed of 
opaque and clear stripes 125 fj,m wide, (b) Blurred image 
obtained conventionally using full illumination and taken in 
a single 50-/is frame; Rayleigh diffraction bound is 1.86 mm. 
(c) Sub-Rayleigh image using focused illumination and A*' — 
23. Details are obscured by a few pixels with very high event 
counts, (d) 3-D intensity profile of (c) with the stripes clearly 
revealed by clipping pixels with very high event counts at the 
Umit of 800. 

To implement sub-Rayleigh imaging in Fig. [TJa) we 
manually scanned the focused beam in a random pat- 
tern, making sure that there was coverage for the entire 
area of interest. At each scan location, we recorded over 
8000 measurement frames for image averaging that took 
less than 1 second to accomplish. The incident power was 
adjusted to have an average peak photocount iVpcak =14 



per integration time (for one pixel). For each measure- 
ment frame, each pixel with exactly N photocounts (after 
dark-count subtraction, measured separately) was tagged 
as having an A'^-photon event. All other pixels were then 
tagged for zero A^-photon events. The measurement pro- 
cess was then repeated at a different scan location un- 
til the object of interest was fully scanned. Figure [2][c) 
shows the resultant image for A^ = 23, revealing the three 
stripes that were lost under conventional illumination in 
Fig.ElJb). The color scale of Fig.[2ljc) has a large range to 
accommodate several pixels with very high event counts 
and therefore image details (with lower event frequencies) 
are obscured. Figure [2][d) shows the 3-D intensity profile 
of the same A^ = 23 image of Fig. [SJc), except that we 
cap the event occurrence at 800 for pixels with over 800 
events in order to make the lower-count pixels more visi- 
ble, thus revealing the three stripes very clearly. Accord- 
ing to theory, the expected enhancement of vW should 
yield a sub-Rayleigh resolution of '^1.86/\/23 = 0.4 mm 
is qualitatively borne out by our results. Note that the 
sub-Rayleigh resolution still exceeds the 106-/xm-limit set 
(after magnification) by the focused illumination at the 
object, as expected. 

We chose A^ = 23 to be substantially larger than iVpcak 
=14 to avoid the "donut-hole" problem. To illustrate 
this issue, Fig. [3] shows images of a point source obtained 
under various measurement conditions. The aperture di- 
ameter of the imaging optics in Fig. [IJa) was set to 3 mm 
with the same overall image magnification of 5.3 so that 
the Rayleigh bound at the image plane (SPAD array) 
was 1.2 mm. We removed the USAF resolution target so 
that the 20-/im-radius focused spot at the object plane 
served as the point source, and no scanning was neces- 
sary for these images. We took ^-^32,000 measurement 
frames, recorded the photocounts at each pixel for each 
frame, binned them accordingly after subtracting dark 
counts, and processed the data. Figure ^a) shows the 
cross section of the Rayleigh-bound image of the point 
source through the 3-mm-diameter aperture obtained by 
including all photocounts to yield an intensity profile av- 
eraged over the ^32000 frames. We measured Apoak ~ 15 
and, as an indicator for the image size, we obtained the 
full-width at half-maximum (FWHM) of ^Imm. 

Figures [3Ib)-(d) are cross-sectional profiles obtained 
by selecting exactly A^-photocounts for A^ = 9, 15, and 
23, that are respectively smaller than, equal to, and 
greater than iVpeak- For A^ < iVpeak in (b), the center por- 
tion of the point-source image usually received more than 
the threshold level A^ and therefore had few exactly A^- 
photon events. On the other hand, the photocounts away 
from the center decrease from Apeak until the photocount 
average matches the threshold A^ where it shows a peak, 
and hence the image has a "donut-hole" shape [23|. For 
N w Apeak in (c), the profile is single-peaked and looks 
slightly narrower and steeper than the Rayleigh-bound 
image in Fig. EJa) which is much larger than the mag- 




of ^VN is in good agreement with theory. The sub- 
Rayleigh technique may find apphcations in which active 
illumination can be implemented to obtain a higher res- 
olution image at the expense of longer acquisition times. 
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FIG. 3: Image cross section of a point source with a modified 
setup of Fig. Wis,) (3- mm aperture) obtained by collecting: (a) 
all photocounts as in conventional imaging with iVpeak ~ 15; 
(b) exactly A^ =9 photocounts and showing a hole where a 
peak should be, (c) exactly A^ =15 photocounts with a peak 
slightly sharper than in (a); and (d) exactly A'^ =23 photo- 
counts with a sub-Rayleigh peak that is sharper than (c). 



nified size of the focused beam of 125 /im (FWHM) at 
the image plane. When TV :^ Npeak, as required for sub- 
Rayleigh imaging, we observe in (d) that the profile for 
A^ = 23 shows a much narrower peak with a FWHM 
width of ^0.4 mm that is smaller than the 1-mm width 
(FWHM) of the Rayleigh-bound point source image in 
Fig. [31a). We also note that the A^ = 23 event frequency 
is much lower than for N less than or equal to A^poak, 
because such a large A^ did not happen very often. Our 
sub-Rayleigh imaging has characteristics that are simi- 
lar to A^-photon interferometry, in which sub- wavelength 
(X/N) interference patterns were obtained with coherent- 
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state input and postselective A^-photon detection 

Our point-source images suggest the following physical 
origin for the formation of the Rayleigh diffraction bound 
and the basis for our sub-Rayleigh imaging technique. 
The donut holes that form for A^ < iVpeak are responsi- 
ble for the diffractive spread of the imaging system. By 
removing the low-A photocount events and those com- 
parable to A^peak, one capturcs the much narrower profile 
for A^ ^ A^peak at the expense of longer acquisition times 
due to less frequent occurrences for high values of A^. We 
should also achieve similar results if we relax the mea- 
surement requirement from exactly A^ for A" ^ Apeak to 
a sum of all A^ » A^pcak, whose resolution should be dom- 
inated by its lowest A^ term due to its larger sub-Rayleigh 
width (oc I/vW) and higher rate of occurrence. 

In conclusion, we have demonstrated sub-Rayleigh 
imaging resolution using a classical light source, tight fo- 
cusing on the object, and A^-photon photodetection. The 
sub-Rayleigh technique removes the low-A^ image com- 
ponents that contribute to the diffractive spread of the 
imaging system. The measured resolution enhancement 
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